Topological winding properties of spin edge states in Kane-Mele graphene model 



Zhigang Wang/ Ningning Hao/ and Ping Zhang^'^'Q 

' LCP, Institute of Applied Physics and Computational Mathematics, 
P.O. Box 8009, Beipng 100088, People's Republic of China 
'''institute of Physics, The Chinese Academy of Sciences, Beijing 100080, People's Republic of China 
Center for Applied Physics and Technology, Peking University, Beijing 100871, People's Republic of China 

We study the spin edge states in the quantum spin-Hall (QSH) efTect on a single- atomic layer 
graphene ribbon system with both intrinsic and Rashba spin-orbit couplings. The Harper equation 
for solving the energies of the spin edge states is derived. The results show that in the QSH phase, 
there are always two pairs of gapless spin-filtered edge states in the bulk energy gap, corresponding 
to two pairs of zero points of the Bloch function on the complex-energy Riemann surface (RS). The 
topological aspect of the QSH phase can be distinguished by the difference of the winding numbers of 
the spin edge states with different polarized directions cross the holes of the RS, which is equivalent 
to the Z2 topological invariance proposed by Kane and Mele [Phys. Rev. Lett. 95, 146802 (2005)]. 
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Since the discovery of the integer quantum Hall effect 
(IQHE) about three decades ago ^|], topological invari- 
ants, which classify the electronic states, have been well 
accepted as a powerful tool for understanding the quan- 
tum many-body phases which have bulk energy gaps ■ 
Historically, the pioneer work was performed by Thou- 
less, Kohmoto, Nightingale, and den Nijs (TKNN) 
who recognized that the IQHE can be understood in 
terms of topological invariants known as Chern num- 
bers 3, which are integrals of the fc-space Berry cur- 
vatures of the bulk states over the magnetic Brillouin 
zone. While IQHE finds its elegant connection through 
the adiabatic curvature with bulk topological invariants, 
Halperin Q first stressed that the existence of the sam- 
ple edges, which produces the current-carrying localized 
edge states in the Landau energy gap, is essential in the 
Laughlin's gauge invariance argument Hatsugai fur- 
ther developed a topological theory of the edge states 
, in which topological invariants are the winding num- 
bers of the edge states on the complex-energy Riemann 
surface (RS). 

Very recently, another striking topological quantum 
phenomenon, i.e., the quantum spin Hall effect (QSHE), 
was identified [1, Q after long-distance efforts on metal- 
lic and conventionally insulating spin Hall effects 
[1^ , and soon after has been attracting extensive current 
interest [11 [11 [11 [M [B [11 H S [111, [12, [11, [11, 
[H, [11] due to its basic physics and its potential applica- 
tion in dissipationless spintronics. Unlike charge IQHE, 
whose presence fundamentally rely on the breaking of the 
time-reversal (T-) symmetry via external magnetic field 
[l| or intrinsic magnetic gauge flux [13, Ull, QSHE does 
not violate the T-symmetry, which implies the absence of 
the non-zero Chern invariants in QSHE insulators. Then, 
it turns out that a Z2-valued topological invariant could 
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be associated with QSHE This Z2 topology is, as one 
selective choice, characterized by whether the number of 
Kramers doublet localized at the edges in a strip geome- 
try is even (nonzero) or odd. If even, the insulating phase 
is an ordinary Bloch insulator; otherwise, the insulating 
phase is a QSHE insulator. 

To study QSHE and Z2 topological order, Kane and 
Mele introduced the model of graphe ne p , which con- 
sists of two copies of Haldane's model [27'|, one for spin- 
up electrons along some axis and one for spin-down elec- 
trons. T-symmetry can be maintained if the intrinsic 
IQHE magnetic fields are opposite for the two spin com- 
ponents. To make the model more physical in that real- 
istic mixing of the two spin components should emerges, 
a Rashba spin-orbit coupling term is further included 
Q , which now becomes the well-known Kane-Mele (KM) 
model. The KM model is the simplest possible model in 
that it has four spin-split bands, which is the minimum 
number required for the nontrivial phase to exist [26j . 
For this reason, the KM model has received the most at- 
tention in the other studies [l^, [H [H, ^2.9] , in which the 
main focus is on the boundary phase twist and disorder 
effect on QSHE. 

In this paper, we give a topological study of the spin 
edge states and its relation with QSHE by using the KM 
model. This study closely parallels with Hatsugai's topo- 
logical theory of edge-state IQHE in that we are seek- 
ing the winding numbers of the spin edge states on the 
complex-energy RS. We show that the spin-edge-state en- 
ergy loops cross the holes of the RS, generating winding 
numbers /| and /| for spin-up and spin-down electrons, 
respectively. The quantized charge (spin) Hall conduc- 
tance could be expressed as a summation (difference) of 
the spin-up and spin-down winding numbers. Thus we 
propose an edge-state topological invariant Is=I]—Ii to 
distinguish a quantum spin-Hall (QSH) insulator from an 
ordinary insulator. If Is is zero, the insulating phase is 
an ordinary insulator; otherwise, the insulating phase is a 
QSH insulator. We stress that this classification between 
topological and ordinary insulating phases survives but 
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FIG. 1: (Color online) Structure of a graphene ribbon with 
zigzag edges, consisting of sublattices A and B. The width of 
the graphene ribbon is A'^. Every unit cell contains A'' numbers 
of A and B sublattices. 



is strongly modified by the weak mixing of the two spin 
components, which reflects the fact that although the ex- 
act quantization of the spin Hall conductance is destroyed 
by the spin-nonconserved perturbation, the QSH phase 
is still topologically distinct from the ordinary insulating 
phase. 

Now we consider the tight-binding KM model of 
graphene [1, 0|, which generalizes Haldane's model [l^l 
to include spin with T-invariant spin-orbit interactions, 

H = tY^ cl^Cja, + iXsO ^y4aSa/3Ci/3 (1) 

{ij)a {{ij))af3 

+ iXR 4a (s X d^j) ^ Cj/3 + A„ ^ ^^cla'^^a■ 

{ij)a i 

The first term is the nearest neighbor hopping term 
on the graphene (honeycomb) lattice, while the sec- 
ond term is the mirror symmetric spin-orbit interac- 
tion which involves spin-dependent second neighbor hop- 
ping. Vij — {2/ \/S){di xd2)z=±l, where i andj are next- 
nearest neighbors, and di and d2 are unit vectors along 
the two bonds that connect i to j. is a Pauli matrix 
describing the electron's spin. The third term, which will 
arise due to a perpendicular electric field or interaction 
with a substrate, is a nearest-neighbor Rashba coupling 
term, which explicitly violates the z—^ ~z mirror symme- 
try. The fourth term is a staggered sublatticc potential 
(^i=±I), which violates the symmetry under twofold ro- 
tations in the plane. 

For the bulk graphene, this Hamiltonian can be written 
in the momentum space. For each k the Bloch wave 
function is a four-component eigenvector |u(k)) of the 
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with X=i sin(k-ai) 

Z=-2AsoE»sin(k-b,) 
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Y=tJ2^cos{k■ai), 
and 

Tifgikai^^ Here 
ai=(^-^, -i^ a, a2=(0, l)a, and a3 = (^^,-i^a 

represent the vectors from sublattice site A (red cir- 
cles in Fig. 1) to its three nearest sublattice sites B 
(green circles in Fig. 1), respectively, and bi=a2— as, 
b2=a3— ai, and b3=ai— a2 represent the vectors be- 
tween the nearest same sublattice sites. No doubt, the 
Hamiltonian ^ can also be written in terms of the 
SO (5) Clifford algebra, which is a very helpful technique 
to deduce Z2-valued topological invariants 30]. 

To explore the edge topological invariant characteriz- 
ing the QSH phase, we now turn to study the graphene 
ribbon with zigzag edges (see Fig. 1), which is periodic 
in the x direction while it has two edges in the y direc- 
tion. In the following, we replace index i with (mns) to 
denote the lattice sites, where (mn) label the unit cells 
and s label the sublattice A and B in this cell. Since the 
graphene ribbon is periodic in the x direction, we can use 
a momentum representation of the electron operator 



ikX,. 



'7ris,a(fc), 



(3) 



where (^(mns) i ^(mns) ) represents the coordinate of the 
site s in the unit cell (mn) , and k is the momentum along 
the X direction. For simplification, let us first consider 
the spin-conserved case, i.e., the Rashba term in Hamil- 
tonian ([T|) is set to be zero. Also, we let the staggered 
sublattice potential vanishes and only consider the first 
two terms in the Hamiltonian ([T]) . Now let us consider the 
one-particle state \i^{k))=J2n,s,a'^risa{k)jl^^jO) . In- 
serting it into the Schrodinger equation H\ip)=e\'ip), one 
can get the following eigenvalue equations for A- and B- 
sublattice sites: 

(e +P2)V'«/lT = PiVAiBT + #(n-l)BT 

+ P3 bl^(n+l)A^ + 1|J{n-l)A^] , (4) 
(e - P2)tlJnAl = Pl^'nBi + t^'(„-l)i3i 

- P3 [*(„+l)Ai + *(„-l)Ai] , (5) 
(e - P2)i^nB] = Vl^nA] + ti^(n+l)A] 

- P3 [V'(n+1)BT + '^{n-l)B]] , (6) 

+ P3 + V'(n-l)i3i] , (7) 

where pi=2t cos (^^ka^ , p2=2Aso sin (\/3fca) , and p^ = 
2Aso sin (^^kaj . Eliminating the i3(A)-sublattice sites, 
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FIG. 2: (Color online) (a) Energy spectrum for the graphene 
ribbon with zigzag edges. The spin-orbit coupling parameter 
in Eq. ((T| is chosen as Aso=0.03t and the last two terms are 
neglected. The bands crossing the gap are spin filtered edge 
states, which are depicted in (b) and (c) by plotting their 
spin (sz) distribution on the lattice sites with A;=0.997r and 
fc=1.0l7r, respectively. To be more clear, the energies used in 
the right panels are denoted by ei and £2 for lower and upper 
spin-degenerate edge states, respectively. 

we obtain the difference equation for A(i?)-sublattice 
sites 

fltpn = h + + /a [V'n+l + ^n-l] , (8) 

where fi=e'^ -pl-pl-t^ - 2pl, f2=pl, h=Pit - 2p2Ps. 
and ij^nAai^nBa) was replaced by i/'n- Eq. ([5]) is the so- 
called Harper equation [3l|. Note that since the spin Sz 
is conserved, the spin-up and spin-down electrons satisfy 
the same Harper equation ([5]). 

Now we represent Eq. ([8]) in the transfer matrix form. 
For this purpose we introduce a new wave function 
which is a linear transformation of the original wave func- 
tion ■0n, 

= ± i^'/'n-l + (9) 

where t± = {^b'^ + A{2 + d) ± 6)/2 with h^-h/f2 and 
d=/i//2. The new wave function Lpn can be written in 
the following transfer matrix form 

( U M(e) ( \ , 

with 

Mie) = ) . (10) 

Although ifn is a transformation of the original wave 
function '0n, here we also can take ipo and (^at as the 



wave functions of two edges, which are linked by a re- 
duced transfer matrix as follows: 

where 

mw = [mw]".(J|; *^). (12) 

All kinds of solutions of Eq. p2)) are obtained by different 
choices of ipo and ipi . Now we study the energy spectrum 
of the graphene ribbon with special attention to the edge 
states. The general boundary condition is 

Vw+i = (fo = 0- (13) 

With Eqs. (HI]), ^ and the boundary condition P^ . 
one can easily get that the solutions satisfy 

M2i(e) = 0, (14) 

and 

(PN = Mii{e)ipi. (15) 

If we use a usual normalized wave function, the state is 
localized at the edges as 

J |Afii(e)| ^ 1 localized at n « 1 (down edge), , . 
\ |Mii(e)| > 1 localized at n fa N (up edge). ^ ' 

Because the analytical derivation is very difficult, we 
now turn to start a numerical calculation from Eq. ([1]) 
with Xji and setting to be zero. Figure 2(a) shows the 
energy bands for the graphene with zigzag edges with the 
intrinsic spin-orbit coupling strength Aso=0-03i, same as 
that in Ref . |8i]. From Fig. 2(a) one can see that there are 
two branches of spin-degenerate dispersed energy bands 
with two pairs of spin states lying in the energy gap. 
These two pairs of spin states cross at the T-invariant 
point fc=7r and are localized at the edges of the graphene 
ribbon. To show the spin localization features of these 
gapless spin states, we plot in Figs. 2(b) and 2(c) the 
spin (z-component) distribution along the graphene strip 
for the two pairs of edge states at fc=0.997r and I.OItt, re- 
spectively. At the left side of the crossing point k=TT [Fig. 
2(b)], the spin-up (-down) edge state with lower energy ei 
is localized at the down (up) edge of the graphene ribbon, 
while the spin-up (-down) edge state with upper energy 
£2 is localized at the up (down) edge of the graphene rib- 
bon. On the other hand, at the right side of the crossing 
point fc^TT [Fig. 2(c)], the spin-up (-down) edge state 
with lower energy ei is localized at the up (down) edge 
of the graphene ribbon, while the spin-up (-down) edge 
state with upper energy £2 is localized at the down (up) 
edge of the graphene ribbon. Thus, if we trace the spin- 
up flow of the edge states when varying the momentum 
k in one period, it can be found that the electrons with 
different spin propagate in opposite directions and that 
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the electronic state of graphene is a QSH state. The spin 
Hall conductivity has been shown as <Js—'^ in units of 
e/47r by different methods. 

Before introducing the winding number of the spin 
edge states on the complex energy RS, let us try to inter- 
pret the above two numerical characteristics. The first 
one is that there are always two energy-degenerate edge 
states appearing localized at opposite edges. To explain 
this issue, we suppose the spin edge states to be expo- 
nentially localized on the boundary with the following 
ansatz [33i. i34i] 



(17) 



where A is a complex number. Inserting Eq. (|17p into 
the Harper equation ([5]), one can easily get the complex 
number A satisfying the following equation 



(A + A-i) -6(A + A-i) - (d + 2) 
That means 



0. 



A + A" 



±t±. 



(18) 



(19) 



It is trivial to see from Eqs. and (|19p that if A is a 
solution, then A~^ is also to be. That means if there is 
an edge state localized near one boundary, at the same 
time there is another edge state localized near opposite 
boundary. Moreover, note that the conservation of the 
spin Sz, i.e., the spin-up and -down electrons satisfy the 
same Harper equation ([8]), one easily find that if A is a 
solution describing an edge state with spin-up localized 
near one boundary, then A~^ is another solution describ- 
ing an edge state with spin-down localized near opposite 
boundary. That is the intrinsic reason for the second fea- 
ture that there are always two edge states with opposite 
spins localized near opposite edges. 

Now we show that similar to charge IQHE, QSHE of 
the graphene ribbon is associated with winding number 
of the spin edge states on a complex energy RS. Within 
the topological edge theory this number is, as well as 
the Chern number, a well-defined topological quantity. 
In order to well understand the winding numbers of the 
spin edge states in QSHE, let us first consider that of 
the edge states in charge IQHE. First we ignore the open 
condition and consider the bulk Bloch function at sites 
with y-coordinate of Ly. For Bloch function, Lp\!f' and 

tp^^^ compose an eigenvector of M with the eigenvalue p. 



M(6 



(b) \ / {b) 

(b) - py<^)\ (b) 



(20) 



In order to discuss the wave function of the edge state, we 
extend the energy to a complex energy. In the following, 
we use a complex variable z instead of real energy e. From 
Eq. (HOI) we get 



p(z) = - A(z)-yA^ 



(21) 



and 



Mn{z) + M22{z)-u; 

-Mil (Zj + M22\Z) + UJ 

where A(z)=Tr[M(z)] and cj=^A2(z) - 4. Clearly, 



detM(e) = 1 



(23) 



since detM(e) = l. From Eq. ((22)) . one can find that 
the analytic structure of the wave function is deter- 
mined by the algebraic function ijj=yj^?{z) - 4. The 
RS of uj—yj /S?{z) — 4 on the complex energy plane can 
be built by the conglutination between different ana- 
lytic brunches. The close complex energy plane can be 
obtained from the open complex energy plane through 
spheral pole mapping. Now let us discuss the analytic 
structure oiuj — a/ A2(2) — 4 on the open complex plane. 
If the system has q energy bands, i.e., 



e e [Ai, A2] , [A2J-1, Asj] , [\2q-1A2q] , (24) 

where Xj denote energies of the band edges and Ai < Xj , 
i < j, then lu can be factorized as 



UJ 



\ 



2q 



A,). 



(25) 



In the present graphene system, there are two energy 
bands, so q=2. The two single- valued analytic branches 
are defined on the plane with two secants down or up the 
bank of which we choose corresponding complex angles, 
which is shown in Fig. 3. 

In order to ensure w(/ij) > {j=l,...,g=q — 1) [Here, 
/ij is the energy of the edge state in the gap [A2j, A2j+i]], 
we divide the two single- value analytic branches based on 
the parity of j. Let us consider the case that the energy 
z lies in the band [A2j-i,A2j] [see Fig. 3(a)]. The left 
side of this energy band is the (j — l)th gap, and the 
right side is the jth gap. When the energy z moves from 
the jth band to the {j — l)th gap (the jth gap) along an 
arbitrary path ci (02), only the singularities A2J-1 and 
X2j have contributions to the variance of the principal 
value of argument of w. On the up bank of the secant, 
we distinguish two branches i?+ and i?^ as the following: 

For j is even, if we set arg(z — A2j-i)=0 and arg(z — 
A2j)=7'', which corresponds to a;(/ij_i) > (a;(/ij) < 
0) when z moves along ci (C2), then the branch _R+ is 
defined as a complex plane with q secants. If we set 
arg(z— A2j-i)=27r and arg(z— A2j)=7r, which corresponds 
to Lo{^ij-i) < (w(/ij) > 0) when z moves along ci (02), 
then the branch R~ is defined as a complex plane with q 
secants. For j is odd, the definitions of i?+ and i?^ are 
reverse to those for j is even. So, if z lies in the jth gap 
from the lowest one on the real axis, 
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FIG. 3: (Color online) (a) The open complex energy plane are 
mapped to the close complex energy plane through spheral 
pole projection, (b) Two sheets with two cuts which corre- 
spond to the energy bands of the graphene nanoribbons. The 
RS of the Bloch function is obtained by gluing the two spheres 
along the arrows near the cuts. 



a {-ly cj > 0, z (real) on i?" (a = +, -), (26) 
and at energies /ij of the edge states we have 

u;(m,) - a {-ly |Afn - M22 , (27) 

where /ij e i?", 0;=+,—. In addition, one can easily 
obtain 

m\-~,V'''''"^'^ , (28) 

1 > 2 for J even ' ^ ' 

where the energy e (on i?") is in the jth gap. 

When the branches R^and R~ on the open com- 
plex energy plane are mapped to the close complex en- 
ergy plane through spheral-pole-projection, one can get 
two single-value analytic spherical surfaces. The RS is 
obtained by gluing the two spherical surfaces at these 
branch cuts. Make sure the ± banks face the =F banks 
of other sphere, (see Figs. 3(b)). Note that there are 
two real axes after gluing. In our model, the genus of the 
RS is g=\ for spin- up (or spin-down) electrons, which 
is the number of energy gaps. The wave function is de- 
fined on the g=\ RS Y.g^i{kx). The branch of the Bloch 
function is specified as w > 0, which we have discussed 
above. With Eqs. ([281), (EH) and ([22]), and using the fact 
when ipLy-i{fJ.j) = for fij g i?" and ipL^^i{fij) ^ for 
fj-j € one can obtain that when the zero point is on 

the upper sheet of RS (i?^), the edge state is localized at 
the down edge; when the zero point is on the lower sheet 
of RS (i?~), the edge state is localized at the up edge. 

In Fig. 4, on the RS T,g=i{kx) of our system the en- 
ergy gap corresponds to the loop around the hole of the 
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FIG. 4: (Color online) RS of the Bloch function corresponding 
to Fig. 2(a). To clearly see the degenerate edge states with 
different spin, we schematically separate their traces with dif- 
ferent colors. The orange and purple curves correspond to 
spin-up and spin-down channels, respectively. 

'^g=i{kx) and the energy bands correspond to the closed 
paths vertical to the energy gap loop on the J^g^i{kx). 
The Bloch function is defined on this surface. For the 
fixed kx there is always g=l zero point at the down- 
edge-state energy /i^. Since there are two real axes on 
the Tig^i^kx), correspondingly, there is g=l zero point 
at the up-edge-state energy . 

The above considerations are for the fixed k^. Now, let 
us consider a family of the RSs "Eg^i^kx) parametrized 
by kx changing in one of its periods. 'Sg^i{kx) can be 
modified by this change. However, all the RSs 'Sg=i{kx) 
with different k^ are topologically equivalent if there are 
stable energy gaps in the 2D energy spectrum. 

On the genus g=l RS, the first homotopy group is 
generated by 2g=2 generators, and f3i, i=l. The in- 
tersection number of the curves (including directions) 0] 
is given by (see Fig. 5) 

I{a„Pj) = S,j. (29) 

Any curves on the RS are spanned homotopically by 
and Pi. When ^j{kx) moves p times around the jth hole 
with some integer t, one has 

C(/i,)«/3f, (30) 

which means 

J(a„C(/ij)) = %. (31) 

When the Fermi energy of the 2D system lies in the 
iih energy gap, the charge Hall conductance is given by 
the winding number of the edge state, which is given by 
the number of intersections I{aj,C{fij)) (= I{C{fj.j))) 
between the canonical loop aj on the RS and the trace 
of fij, i.e., 

= -^/(C(m,)). (32) 
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Similarly, the above expression can be obtained by the 
Byers-Yang [sl] and Laughlin-Halperin Q gauge ar- 
guments. The system with the periodic boundary con- 
dition in the x direction and open condition in the y 
direction can be considered as a cylinder. By Laughlin 
gauge invariance argument, the vector potential A has to 
have the form A=njj^{n is an integer). When the flux 
$ threading the cylinder is adiabatically turn on from 
$(0)=0 to $(r)= hc/e with A^—hc/e a flux quantum, 
AA—j^ accords with the gauge argument, therefore, 
Ay4(A$) maps the system back to itself. Basing on the 
single-electron assumption. When the Fermi energy lies 
in the jth energy gap, there are /(C(/Zj)) states (elec- 
trons) transferring from the down edge (?/=l) to the up 
edge {y=Ly — 1) in net. The energy change during the 
adiabatic process is AE=I{C{fij)){—e)V, where is a 
voltage in the y direction. This gives the charge Hall 
current I„ as follows 



Then we get an expression for a°y^° as Eq. p2|) . 

The above analysis is for the winding number of the 
edge states in IQHE and it can be easily generalized to 
the winding number of the spin edge states in QSHE 
when the spin degree of freedom is considered. In the 
latter case, the winding number of the spin edge state 
is given by the number of intersections Is{aj , C'{fij)) 
(= Is{C{iJLj)), 5 =1, J, for spin- up and -down, respec- 
tively) between the canonical loop Uj on the RS and the 
trace of /ij . From Fig. 4 one can observe that fj, moves 
one time across the hole 0=1), which means C(/x) « (3 
and |/s(C(/i))|=l. Considering the winding direction 
(see Fig. 5), one obtain /|(C(//))=1 for spin-up elec- 
trons while /|(C(/i))=— 1 for spin-down ones. Thus, the 
charge Hall conductance is given by the summation over 
two spin channels, 

^^ = aT+^i)X^4^-0, (33) 

while the spin Hall conductance is given by the difference 
between them. 

When the flnite on-site energy is considered while 
the Rashba coupling is kept as zero, the two degenerate 
energy bands will split. For the bulk graphene H(k) 
one can easily obtain the four dispersion energy bands, 

Ck = ±^X^ + Y^ + {Z±Kf, (35) 

between which there is an energy gap with magnitude 
A=2|A„ - 3V3Asol- We plot in Fig. 6 the energy 
dispersion with different A^ and Aso- The condition 
\v 7^ 3\/3Aso provides a flnite bulk energy gap in the 




AAA 
I(A,B)=+1 I(A,B)=-1 I(A,B)=+1 



FIG. 5: Intersection number I[A,B) of two curves A and B. 
Each intersection point contributes by 4-1 or —1 according to 
the direction. 




FIG. 6: (Color online) Energy spectra of the bulk graphene 
(a) in the QSH phase with \v=Q.lt and Aso=0.05t and (b) 
in the ordinary insulator phase with \^=QAt and Aso=0.05t. 
There are bulk gaps appearing in both two different phases. 



graphene. However, when the Fermi energy lies in the 
gap, the graphene's phases are different with different 
topology: for A^, > 3V3Aso: the system is an ordinary 
insulator and for A^j < 3-\/3Asoj the system is a topolog- 
ical insulator with QSH phase. Kane and Mele Q have 
proposed a Z2 index describing the QSH phase, which 
is determined by counting the number of pairs complex 
zeros of P(k)=Pf[(iti(k)|e|Mj (k))] with |ui(k)) the wave 
functions corresponding to the bulk Hamiltonian ([2|) and 
time reversal operator. To understand these different 
phases, Kane and Mele also studied the edge states of the 
graphene ribbon with zigzag edges. However, the wind- 
ing properties of the edge states in the KM model has 
not been previously considered, which just is the special 
focus of our present study. 

Considering nonzero A^;, now let us investigate the 
topological winding numbers with zigzag edges. In this 
case, the Harper equation for site A turns to have the 
same form as Eq. ([8]) , provided that the P2 in the coef- 
flcients /i and /a is replaced by p2 ± \v for spin-up and 
-down channels, respectively. Similarly, one can intro- 
duce a new wave function, as a linear transformation of 
the original wave functions ([9]) , satisfying the matrix form 
with the transform matrix M{e) ITU)) . The p2 in the coef- 
flcients &, d, and t± is now also replaced by p2±A,j for dif- 
ferent spins. In evidence, the edge-state energies become 
non-degenerate for the finite on-site energy A,j . Figure 7 
shows the energy spectrum of the "zigzag" graphene rib- 
bon in a QSH phase with Ai,=0.1i and Aso=0.05i, which 
satisfy A^, < 3\/3Ago. From Fig. 7, one can clearly see 
that the gapless edge states are localized near the system 
boundary with different spins. The corresponding com- 



1 2 



FIG. 7: (Color online) Energy spectrum of the "zigzag" 
graphene ribbon in the QSH phase with Xv=0.1t and 
Aso=0.05t. The red and blue lines represent the edge states 
localized at the down and up edges of the system, respec- 
tively. And the circle and triangle label the up and down 
spins, respectively. 




FIG. 8: (Color online) Energy spectrum of the "zigzag" 
graphene ribbon in the ordinary insulator phase with Ai,=0.4t 
and Aso=0.05t. The red and blue lines represent the edge 
states localized at the down and up edges of the system, re- 
spectively. And the circle and triangle label the up and down 
spins, respectively. 



plex energy RS of the Bloch wave functions is topologi- 
cally equivalent to Fig. 3. And no doubt this QSH phase 
can be described by the difference of the winding numbers 
of the edge states with different spins, i.e., /s=/| — /|=2, 
in units of e/Air. On the other hand, one can obtain the 
charge Hall conductivity CTc is zero since As a 

comparison, we draw in Fig. 8 the energy spectrum of an 
ordinary insulator phase with Ai,=0.4t and Aso=0.05i. 
One can find that in this case (see Fig. 9), there are 
no gapless edge states connecting the two energy bands. 
The traces of the edge states in the corresponding com- 
plex energy RS are sunk in the bulk bands, i.e., /|=/j^=0, 
which means that the system in this case is an ordinary 
insulator without QSHE. We thus conclude that while 
the previouslystudied topological winding number index 
Ic 0, m, [11] has been verified to describe the 

IQH insulator, our present studied topological winding 
number index Is (=/| ~ ) can be used to characterize 
the QSH insulator. 

Finally let us turn to consider the most complex case of 
Sz being nonconserved, i.e., the Rashba coupling A^, as 



FIG. 9: The RS of the Bloch function corresponding to Fig. 
8, in which the edge states are sunk in the energy bands and 
can not form a loop around the hole. 

well as the on-site energy A^ , is finite in the KM graphene 
model (m . After a straightforward derivation, we obtain 
the following couple eigenvalue equations for A_r ^ 0: 

P2±)^nA± — Pl'^nB± +^l^(n-l)_B± 

±P3 [*(ri+l)A± +*(n-l)A±] 

(eTP2±)*riB± =Pl'i?nA± +tl1'(„+l)A± 

TP3 [* (n+l)B± + *(n-l)B±J 

- ip±'i/nAT + «Afl*(„+i)A=F, (37) 

where p2±=P2 ± A„, p±=2Ai^ cos(^fca± -1), and "±" in 
'^nA{B)± label the spin channels. Although the Harper 
equation and the corresponding transform matrix are 
very difficult to derive in this case, we can also distinguish 
the QSH phase from the ordinary insulator phase by the 
topological winding index introduced in the present pa- 
per directly from the energy spectrum. As an example, 
we reinvestigate the QSH phase and the ordinary insula- 
tor phase with the Rashba coupling being a finite value 
Xu=0.1t. The other parameters are chosen be the same 
as thosed used in Fig. 7. The calculated energy spec- 
trum (not shown here) turns to be essentially the same 
as that depicted in Fig. 7. As a result, the corresponding 
RS of the Bloch wave functions is topologically equiva- 
lent to Fig. 3. This means that in spite of finite Rashba 
term, the system is still in the QSH phase with the topo- 
logical index Is—2. Therefore, although in this case the 
spins are nonconserved (see Table I) and the exact quan- 
tization of the spin Hall conductance is destroyed by the 
spin-nonconserved perturbation, the QSH phase of the 
graphene system is still topologically distinct from the 
ordinary insulating phase once provided that the Rashba 
term does not change the system's topological proper- 
ties. This conclusion keeps consistent with Ref. Q. If 
the Rashba spin-orbit coupling Xn turns much strong, the 
concept of the "spin edge states" becomes faint and the 
winding number of the spin edge states becomes mean- 
ingless. 

Before ending the present paper, we would like to 
stress that the topological index Ig describing the QSH 
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TABLE L Comparison of the spin probability distribution 
((s^)) of the lowest-energy edge states on the lattice sites in 
the cases Xr—Q and \u=0.1t. The other parameters are set 
as Aso=0.05t and A„=0.1f and the momentum k is kept as 
fc=0.997r. 



Case Ah=0 


Case A_R=0.1t 


lattice site index n 




lattice site index n (sz) 


1 


0.9998 


1 0.9635 


2 


0.0001 


2 -0.0274 


3 


0.0001 


3 0.0001 


Total probability 


1.0000 


Total probability 0.9362 




with topological integer Is=I] — Ii, there are /| spin-up 
electrons transferred from one edge to the other when a 
unit magnetic flux is adiabatically through a cylindrical 
system by /| branches of gapless edge states. And at the 
same time there are j/J — /| spin-down electrons to be 
transferred with the opposite direction (see Fig. 10). 

In summary, by theoretically studying within the KM 
Hamiltonian a graphene strip with zigzag edges, we have 
presented an alternative topological index Is to charac- 
terize the QSH phase. The topological index Ig describ- 
ing QSH phases has been defined as the difference of the 
winding numbers of the spin- resolved edge states crossing 
the holes of the complex-energy RS. Based on this topo- 
logical index, we have discussed different phases by mod- 
ulating different parameters in KM model, which agree 
well with the previous studies in terms of the conven- 
tional Z2 topological invariance. 



FIG. 10: (Color online) The Laughlin-Halperin diagram of 
the QSH system. 



phase can also be interpreted in terms of Laughlin and 
Halperin's arguments |5j, 16)]. In fact, in the QSH system 



Acknowledgments 

This work was supported by NSFC under Grants Nos. 
10604010, 10534030, and 60776063, and the National 
Basic Research Program of China (973 Program) under 
Grant No. 2009CB929103. 



[1] K. V. Klitzing, G. Dorda, and M. Peper, Phys. Rev. Lett. 

45, 494 (1980). 
[2] X. G. Wen, Phys. Rev. B 40, 7387 (1989). 
[3] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. 

den Nijs, Phys. Rev. Lett. 49, 405 (1982). 
[4] D. J. Thouless, Topological Quantum Numbers in Nonre- 

latisvistic Physics (World Scientific, Singapore, 1998). 
[5] B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 
[6] R. B. Laughlin, Phys. Rev. B 23, 5632 (1981). 
[7] Y. Hatsugai, Phys. Rev. Lett. 71, 3697 (1993); Phys. 

Rev. B 48, 11851 (1993). 
[8] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(2005). 

[9] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

[10] S. Murakami, N. Nagaosa, and S.C. Zhang, Science 301, 
1348 (2003). 

[11] J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jung- 
wirth, and A. H. MacDonald, Phys. Rev. Lett. 92, 126603 

(2004) . 

[12] S. Murakami, N. Nagaosa, and S.-C. Zhang, Phys. Rev. 

Lett. 93, 156804 (2004). 
[13] M. Onoda and N. Nagaosa, Phys. Rev. Lett. 95, 106601 

(2005) . 

[14] B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 
106802 (2006). 

[15] X.-L. Qi, Y.-S. Wu, and S.-C. Zhang, Phys. Rev. B 74, 
085308 (2006). 

[16] D. N. Sheng, Z. Y. Weng, L. Sheng, and F. D. M. Hal- 
dane, Phys. Rev. Lett. 96, 036808 (2006). 



[17] T. Fukui and Y. Hatsugai, Phys. Rev. B 75, 121403(R) 
(2007). 

[18] L. Fu and C. L. Kane, Phys. Rev. B 74, 195312 (2006). 
[19] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

[20] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 

[21] S. Murakami, Phys. Rev. Lett. 97, 236805 (2006). 

[22] M. Onoda, Y. Avishai, and N. Nagaosa, Phys. Rev. Lett. 

98, 076802 (2007). 
[23] C. Wu, B. A. Bernevig, and S.-C. Zhang, Phys. Rev. Lett. 

96, 106401 (2006). 
[24] C. Xu and J. E. Moore, Phys. Rev. B 73, 045322 (2006). 
[25] B. A. Bernevig, T. L. Hughes, S.-C. Zhang, Science 314, 

1757 (2006). 

[26] J. E. Moore and L. Balents, Phys. Rev. B 75, 121306 (R) 
(2007). 

[27] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 
[28] K. Ohgushi, S. Murakami, and N. Nagaosa, Phys. Rev. 

B 62, R6065 (2000). 
[29] A. M. Essin and J. E. Moore, Phys. Rev. B 76, 165307 

(2007). 

[30] S. Murakami, N. Nagaosa, and S. C. Zhang, Science 301, 
1348 (2003); Phys. Rev. B 69, 235206 (2004). 

[31] P. G. Harper, Proc. Phys. Soc. London Sect. A 68, 874 
(1955). 

[32] N. Byers and C. N. Yang, Phys. Rev. Lett. 7, 46 (1961). 
[33] M. Konig, H. Buhmann, L. W. Molekamp, T. L. Hughes, 

C.-X. Liu, X.-L. Qi, and S.-C. Zhang, arXiv/0801.0901. 
[34] M. Creutz and I. Horvath, Phys. Rev. D 50, 2297 (1994); 

M. Creutz, Rev. Mod. Phys., 73, 119 (2001). 



9 



[35] Z. Wang and P. Zhang, Phys. Rev. B 77, 125119 (2008). Phys. Rev. B 78, 075438 (2008). 

[36] N. Hao, P. Zhang, Z. Wang, W. Zhang, and Y. Wang, 



